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The Weyl-gauge {Aq — 0) QCD Hamiltonian is unitarily transformed to a representation in 
which it is expressed entirely in terms of gauge-invariant quark and gluon fields. In a subspace 
of gauge-invariant states we have constructed that implement the non-Abelian Gauss's law, this 
unitarily transformed Weyl-gauge Hamiltonian can be further transformed and, under appropriate 
circumstances, can be identified with the QCD Hamiltonian in the Coulomb gauge. We demonstrate 
an isomorphism that materially facilitates the application of this Hamiltonian to a variety of physical 
processes, including the evaluation of 5-matrix elements. This isomorphism relates the gauge- 
invariant representation of the Hamiltonian and the required set of gauge- invariant states to a 
Hamiltonian of the same functional form but dependent on ordinary unconstrained Weyl-gauge 
fields operating within a space of "standard" perturbative states. The fact that the gauge-invariant 
chromoelectric field is not hermitian has important implications for the functional form of the 
Hamiltonian finally obtained. When this nonhermiticity is taken into account, the "extra" vertices 
in Christ and Lee's Coulomb-gauge Hamiltonian are natural outgrowths of the formalism. When 
this nonhermiticity is neglected, the Hamiltonian used in the earlier work of Gribov and others 
results. 



I. INTRODUCTION 



In earlier work on QCD in the Weyl gauge {Aq = 0), we have constructed gauge-invariant operator- valued quark 
and gluon fields ;(jy these include the gauge- invariant quark field 

V'Gi(r) =Fc(r)V(r) and Vl,(r) - V^t(r) i/-i(r) , (1) 

where 

Vc{r) = exp(-*53^(r)^) exp (-zgA'^(r)4) , (2) 



yc"i(r)=exp(i5A"5(r)4) exp ( igiV" (r) ^ ) , (3) 

and where the A" designate the Gell-Mann matrices. In these expressions X"{r) = [^A^(r)], so that diX"{r) is the 
i-th component of the longitudinal gauge field, and 3^"(r) is defined as 3^"(r) — [ ^Ajir)]. Ajir), which we refer 
to as the "resolvent field" , is an operator-valued functional of the gauge field, and is represented in Refs. and Q 
as the solution of an integral equation. Constructing a gauge-invariant quark field by attaching Vc(r) to the quark 
field tp represents an extensioii, into the non-Abelian domain, of a method of creating gauge-invariant charged fields 
originated by Dirac for QED; Q and, like Dirac's procedure, this non-Abelian construction is free of path-dependent 
integrals. An explicit demonstration that ^Gi(r) is invariant to non-Abelian gauge transformations has been given 
by implementing gauge transformations with the generator exp{— zJdyCJ°(y)cL;°(y)} where is the non-Abelian 
"Gauss's law operator" 

G'^ = d,nt+ gr'^Aln^ + #t (r) f ^(r) , (4) 
and a;° is a number-valued gauge function. With the use of this generator, under which 

^(r) ^ ^'(r) = exp {-iu;"{r) ^ ) V(r) (5) 

and 

Al{r)^ ^ exp(-z^"(r)^) (^.^(r)^ -f i^,) exp (z^"(r)^ ) , (6) 

it has been shown that Vc (r) also gauge-transforms as 

Vc(r)^Fc(r)exp(ic^"(r)^) and Vc-\r) ^ exp {-iu;''{r) ^) Vc\r) (7) 
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so that ■0Gi(r) remains gauge- invariant. P The resolvent field Aj also has an important role in the gauge-invariant 
gauge field 



AGH(r) = [A^„(r) ^] = Vc{r) [At{r) >^]V,-\r) + i Vc{r) d.V,-\r) , 
which can be shown to be the transverse field Q 



S- ■ - Ml 



A){v) + A'^{v) 



(8) 



(9) 



Eq. |(Sll, as well as the fact that ^Qj(r) and ^"^(x) commute, demonstrate that A'q^^[t) is gauge- invariant — more 
precisely, invariant to "small" gauge transformations. We can also define a gauge-invariant chromoelectric field 



'IIq j. 5\ A natural definition of IIq ^ in this formulation is 



nci.W = [n^H(r) ^] = Vc{r) ^ Vc-\r) n^(r) 



or, equivalently, 



RabTl' where Rab = hM>^"VcX''Vc-^] 



(10) 



(11) 



where 11° is the momentum conjugate to the gauge field Af in the Weyl gauge. With the use of the commutator 

"^^(x), 7^afc(y)] = ^.g/^''«7^a,(y)(5(x - y), (12) 

obtained in Ref. 0, it is easy to verify that nQj(y) commutes with G'^ix) and therefore also is gauge- invariant. 

In this work we will use a representation, which we discuss in Section^ in which the Weyl-gauge QCD Hamiltonian 
is expressed entirely in terms of gauge-invariant fields. Since the gauge-invariant gauge field is transverse, it is of 
interest to relate this gauge-invariant formulation to the Coulomb gauge. We address this question in Section III Bl In 
Section^lwe also show that the Weyl-gauge QCD Hamiltonian in this representation — in which all operator-valued 
fields are gauge-invariant — must be applied to a set of gauge-invariant states that are solutions of the non-Abelian 
Gauss's law. In Section lTTll we address the problem that these states, which solve Gauss's law in QCD, are complicated 
constructions that are difficult to use. We demonstrate an isomorphism in this section between this Hamiltonian, 
which operates on gauge-invariant states, and a corresponding Hamiltonian that is a functional of gauge-dependent 
Weyl-gauge fields and that operates on a set of "standard" perturbative states. Also, in Section lTTll we relate these 
Hamiltonians to those obtained from Coulomb-gauge formulations of QCD. We discuss the implications of our work 
in Section ITVI 



II. RELATION OF THE GAUGE-INVARIANT REPRESENTATION OF THE WEYL GAUGE TO THE 

COULOMB GAUGE. 



The QCD Hamiltonian in the Weyl gauge has been expressed in terms of gauge-invariant operator-valued fields. 
In this work, extensive use has been made of the unitary equivalence of — the "Gauss's law operator" given in 
Eq. which imposes the non-Abelian Gauss's law — to the "pure glue" version of that operator 



as shown by 



where 



(13) 



(14) 



Uc = exp 



exp 



(15) 



This unitary equivalence has been used to establish a new representation — the Af representation in which 
represents the complete Gauss's law operator Q°', and ip represents the gauge- invariant quark field because it commutes 
with 5°. The Af representation is unitarily equivalent to the C representation in which Q"^ and ipa designate the Gauss's 
law operator and the gauge-invariant spinor (quark) field respectively. In the Af representation, Jq (r) = g^^ (r)^^p{r) 
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and jf (r) = gip''' {r)ai^'ip{r) are the gauge-invariant quark color charge and quark color current densities respectively. 
The Weyl-gauge QCD Hamiltonian can be transformed from its familiar C-representation form 

jdv{ lnUv)m^) + \F^,{v)F^^{v)+^Hv)[Pm-za.{d,-zgA'l{v)^)]i,{v)] (16) 
to the M representation, as shown by 

Hg,=Uc-^HUc. (17) 

This similarity transformation leaves the gauge field untransformed, but it transforms the quark field and the negative 
chromoelectric field as shown by j^] 

Z^c~'(x) V'(x)Z^c(x) = Fc-i(x)^(x) (18) 

and 

Uc-'{^)Iin^)i^ci^) = n?(x) - Rbai^)dt^J dyI?''^(x,y)jS(y) . (19) 

The transformed, A/"-representation Hamiltonian H can be expressed entirely in terms of gauge-invariant variables by 
making use of the identities 

RaqRbq = Sab, f'^'^RuaRvb = and d.Rba = - F'"'' RuaPvr whcrc P„ = -^Tr[A''Fc d^Vc^] . 

The QCD Hamiltonian in the M representation, expressed in terms of gauge-invariant fields, is 

^Gi = ! dv 



where 



iH^,Vr)Hgn(r) + \F^u,{^)F'^^^J{^) " (/3m - la^di) V'(r) 

i j dKdy (^J„"(;,)(x)P'^''(x,y)jo^(y) + jo^(y) 2?^" (y, x) Jo"(g,) (x)) - 

i J drdxdyj§(y) V^^ (y, r)d'V'^''{r, x)jo'(x) - J drj7(r)Ag| ,(r) + Hg . (20) 
i^GV.,(r) = a,Agu(r) - 9.^a,^.(r) - 5r^^A^„(r)Aa,,(r) , (21) 



from which it follows that 

FS,,^ir) = Ra,{r)F^^{r). (22) 

Because Hq\ is in the Af representation, ip and ■0^ denote the gauge-invariant quark fields. I?°''(x, y) is the inverse 
Faddeev-Popov operator, which we will discuss in Section III Al and Jq(^q\^{'ic) is the gauge-invariant gluon color charge 
density, defined as 

^0(Gi)W-5r'^4u(r)naH(r). (23) 

T<^t ^ Tib u . n„A nb , 



Although Hq\ is hermitian, Hq ^ is not, because, as can be seen from Eq. ((TT|l . Hg, - — U^Rab, and H^ does not commute 
with Rab- Similarly, Jq^^i-, is not hermitian, and >/o(Gi) ~ 9 f'^'"^^'G\ i^G\ i ■ '^^'^ ^^^^ P^'^^ of the QCD Hamiltonian is 

ffe = dxdy[eg,(x)I?'^^(x,y)jo^(y) + jS(y)P'"^(y,x)gg,(x)] (24) 

where C/q is the gauge-invariant Gauss's law operator 

which consists solely of gauge-invariant fields, every one of which commutes with Q"', the Gauss's law operator in the 
Af representation; ^/q is hermitian because i?., and commute. @. 

Eq. (|2()|l resembles the QCD Hamiltonian in the Coulomb gauge. The only direct interaction between color currents 
jf and the gauge field involve the transverse current only. The other interactions in which quarks participate are 
nonlocal, involve the quark color-charge density Jq , and are mediated by Green's functions that are the non-Abelian 
generalizations of the Abelian d~^. These interactions still involve the longitudinal component of the gauge- invariant 
chromoelectric field, but we will show how this can be eliminated in Section IH Bl 



A. The inverse Faddeev-Popov operator. 

The Faddeev-Popov operator in the gauge-invariant representation of the Weyl gauge is 

di and Di commute because Aq^^ is transverse. The Faddeev-Popov operator has a formal inverse, which can be 
represented as the series 

oo 

P^''(y,x) = ^ /(^f {Tl^iy)6{y - x)) , (26) 

n=0 

where f"^-!'- represents the chain of SU(3) structure constants 

ySh'i _ jaibsi jsia2S2 ^^20383 ^ ^ ^ JS(„_2)a(„_i)S(„_i) jS(„_i)a„?i (27) 

and where repeated superscripted indices are summed from 1^8; for n ~ 1; the chain reduces to /f'''' ~ and 
for n = 0, Jq''''' = —Sbh- ^")(r)jo C'^) ^ special case of a general form T^"^(r)(p'*(r) for an arbitrary (p'*(r) given by 

%(r)^^r) = <^|,)(r) %1 «|,)(r) %1 (• • • «t)(r) %^ (^''(r))))) , (28) 

with 

T^%{rW\r) = ^''(r) and 7^f)(r)^"(r) = Ag,,(r) |i^''(r). (29) 

By expanding I?''''(y,x) and combining terms of the same order in g, it can be observed that, as will be proven in 
Appendix EI 

d-D^^^V''''iy,^)^Sab5{y-^) (30) 

where D"'' = diSah + g/^'^'^A^u and that 

2?^''(y,x) 9.P'-(,)=<5b,%-x), (31) 

where 

(32) 

and 

d^'''>= (^92 s^, - gfib d, Al, ^ (33) 

and the ^ symbol indicates that and di differentiate to the left. In demonstrating Eqs. H3UI) and H31|) . it can be 
helpful to use the expanded form of the n-th order term of the inverse Faddeev-Popov operator series 

P(„)(y,x) -,g/ / / J i^\y^,i,)\'^OuM'))dz{i)J An\z{i)~z{2)\'' 

^^^'^^^^''>d^/--J 4^|z(„-i)) - z(„)r<=""^"^"^^az(„),„ 4n\z{n) x| ^^^^ 



with 
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and 



iT/Wy-x) = ./"•7 (j^^) . (36) 

Integration by parts with respect to the z{i) and the identity f^^-j^ = i^^Y'' demonstrate that 

V'^\y,^)=V'^^^,y). (37) 
It is apparent from Eqs. that 'D^^{y,y.) obeys the integral equation 



which has these equations as an iterative solution. 

Eq. H2ti|l enables us to express the commutator of the gauge-invariant gauge field and the negative gauge-invariant 
chromoelectric field as 

K„(y),^GHW] - {5^b5^A^-y)) + ^P'"'(y,x) D^- (x)) . (39) 
Eq. H39(l and the commutator, obtained in Ref. [3, 



[ng„(x) ,n^„(y)] = |^i?'^''(x,y)/''^^n^,^.(y) - Ap'">(y,x)/"-n^,,(x)| , 



(40) 



are in agreement with those given by Schwinger for the Coulomb gauge, except for some differences in operator order. 
This fact suggests that the gauge-invariant Weyl-gauge field and the Coulomb-gauge field discussed by Schwinger are 
very similar. The differences in operator-order should be expected because, in Ref. [9j, ambiguities in operator order in 
the Coulomb gauge are resolved by symmetrizing noncommuting operator-valued quantities so that Coulomb-gauge 
operators are kept hermitian. In our work in the gauge-invariant formulation of the Weyl gauge, ambiguities in 
operator order do not arise. When, because of a non-symmetric ordering of gauge fields and chromoelectric fields, 
some gauge-invariant operator-valued quantities turn out not to be hermitian, we leave them that way in order to 
avoid ad hoc changes in operator order. 

Eq. H4U|I leads to the commutation rule for the transverse parts of IIq ^ (y) , 

[ngT(x),n^T(y)] =0. (41) 

Eq. (|39(l leads to the commutator of the transverse part of nQ^ (y) and ^(x) (which is transverse) 

[n^b(y)'^Gu(x)] =-*^a6(%-^)'5(x-y); (42) 
Eq. (|39l) can be shown to be consistent with 9^ ^ = because 



d 

dxi 



[n^„(y),^Q.(x)] = -^(SatS.,^Si^-y) + ^V''^{y,^)^^^^^ 



-^SabS^, ( ^<5(x - y) + ^^(x - y) ) - , (43) 



and with DjUQ^-fnO because 

^f(y)N„(y),^Gu(x)] =0 (44) 

trivially. 

The Faddeev-Popov operator has a well-documented importance in non-Abelian gauge theories. Gribov has shown 
that gauge fields that have been gauge-fixed to have a vanishing divergence can differ from each other, [Til and 
that the Faddeev-Popov operator does not have a unique inverse. In that same work, Gribov makes the suggestion 
that the zeros of the Faddeev-Popov operator d^Sac + gf°'^'^A\di might so intensify the interaction between color 
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charges that the effect could account for confinement. Subsequent authors have reiterated this suggestion, 0, 0| 
and connections between the zeros of the Faddeev-Popov operator and color confinement have been discussed by other 
authors as well. [IE EE 113 

Eqs. and are based on a series representation of the operator- valued 'D^'^{y,x); they are obtained by 

combining all terms of equal order in g and noting cancellations within each order. They do not, however, establish 
that I?^''(y,x) is the unique inverse of the Faddeev-Popov operator. Questions about uniqueness can readily be 
formulated about number-valued functions, but are very difficult to address for operator-valued quantities. Eqs. H3U|) 

and establish that 2?^''(y,x) is an operator- valued inverse of (acting on the left) and of d-D'^^-j (acting on 

the right) without addressing the question of its uniqueness. However, although ^ is an operator- valued quantity, 
the SU(2) versions of its constituents — the Weyl-gauge field Af and the resolvent field Af — can be, and often 
have been, represented by number-valued realizations as functions of spatial variables. Such realizations have been 
used extensively to study the topology of gauge fields. 0,0,0] When the integral equation for the resolvent field 
referred to in Sectionals expressed in terms of a number- valued hedgehog representation, it can be transformed into 
a nonlinear differential equation that was shown to have multiple solutions. ^ Moreover, this nonlinear differential 
equation was shown to be very nearly identical in form to the one used by Gribov as a specific illustration of the 
fact that the Faddeev-Popov operator for the transverse SU(2) gauge field does not have a unique inverse. With this 
number-valued realization we were able to establish that the gauge-invariant field, which is transverse, has a Gribov 
ambiguity, Q even though there are no Gribov copies of the gauge-dependent Weyl-gauge field. !l9, 20, 21] 

In the context of the quantized theory — for example, in Hq\ — we will represent 'D^^{y, x) as the operator- valued 
series described in Eqs. H26(l and l|34|l . Since each term in this series has unambiguous and self-consistent commutation 
relations with all other operator-valued quantities, the series representation of I?''''(y,x) is entirely satisfactory for 
determining the commutators of Hq\ with other gauge-invariant operators — and therefore determining their time 
dependence — even though number-valued realizations of the gauge-invariant gauge field lead to nonlinear integral 
equations that do not have unique solutions. 

It may seem surprising that, starting in the Weyl gauge and expressing the QCD Hamiltonian in that gauge in 
terms of gauge-invariant variables can lead to a form of the Hamiltonian that, while never actually having been gauge- 
transformed, has the same dynamical effect as the QCD Hamiltonian in the Coulomb gauge. But a remarkably similar 
state of affairs obtains in QED. When QED is formulated in the temporal gauge, and a unitary transformation is carried 
out that is the Abclian analog of the one that leads to the Hamiltonian described in Eqs. H20l) - (|24|l . the following result 
is obtained: (22. .2,1] The QED Hamiltonian in the temporal gauge, unitarily transformed by exp \i J ■^diAi{r)jQ{r)dr\ 
— the Abelian analog of the transformation Uc described in Eq. H15|l — takes the form 

Hqed ^ fdr [iH,(r)H,(r) + lF,^{r)F,j{r) + i;\r) {fim - iaA) ^{v)] 



dvMv)Al{v) + j dvdv'^-^^^ + H, ■ (45) 

Aj designates the transverse Abelian gauge field — which, in Abelian theories, is also the gauge-invariant field — 
and Hg can be expressed as 



H9 = -\jdv (^a,H,(r)^jo(r) + jo(r)^a.H,(r) 



(46) 



Hg is the Abelian analog of iJg, described in Eq. (|24|l . The Abelian Gauss's law operator, Q — dilii + jo, transforms 
into djli in the representation in which ijj represents the gauge-invariant electron field; and the states that implement 
Gauss's law, which originally are selected by ^(r) ]^'(r)) = 0, are given by 3iHi(r)]$) = in the transformed 
representation (or, as is more appropriate for Abelian gauge theories, by Q^^\y)\^) — and c}iH|^''(r)]$) = 
respectively, where designates the positive- frequency parts of operators). [23, 13 As can be seen, Hqed also 
consists of two parts: the Hamiltonian for QED in the Coulomb gauge; and Hg , which has no effect on the time 
evolution of states that implement Gauss's law, but which "remembers" the fact that Hqed is the transformed 
Weyl-gauge Hamiltonian by preserving the field equations for that gauge. An identical transformation applies to 
covariant-gauge QED, the sole difference being in the form of the Hg produced by the transformation. 

As we can see from Eqs. H20|) . (|24|) . H45|) and 14611 . and as will become even more evident in Eq. H60I) . QCD and QED 
are strikingly similar in the relation between their Hamiltonians in different gauges when these are represented in terms 
of gauge-invariant fields. Nevertheless, there are important differences between QED and QCD in the significance of 
this relationship. One such difference is that, in QED, we may safely use the original untransformed Weyl-gauge or 
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covariant-gauge Hamiltonian in a space of perturbative states when evaluating S'-matrix elements, even though these 
gauge-dependent perturbative states fail to implement Gauss's law. This means that, for perturbative calculations in 
QED, we can safely use the Lagrangian 

C - -iF^.F^'^ - + ^(^7^9^ - m)V^ + L, (47) 

with Lg — — AqG for the Weyl gauge and Cg = —Gd^Af^ + i(l — j)G^ for the covariant gauge, without paying any 
attention to Gauss's law whatsoever. A corresponding practice in Weyl-gauge QCD is the use of the Weyl-gauge 
Hamiltonian i? in a Fock space of perturbative states that are not annihilated by Q'^. There is, however, the following 
important difference between QED and QCD. The use of perturbative states in QED without implementing Gauss's 
law is permissible because, in QED, a unitary equivalence can be established between 9^11^ and 9^11^ + jo, so that 
diHi can be interpreted as diHt + jo in a new representation. |22l In this way, it can be shown that perturbative 
states that implement only dilli{r) « instead of 9^11^ + jo(r) « may be used when evaluating S'- matrix elements in 
QED; the only effect on 5- matrix elements from this substitution consists of changes to the renormalization constants, 
which are unobservable. |25i] But this dispensation to ignore Gauss's law in perturbative calculations has not been 
shown to extend to QCD, because DiHf + jo(r) is unitarily equivalent only to DiHf, but not to 9^11^; and states that 
implement the Gauss's law Dillf w cannot be perturbative states. In particular, the use of Hq\ for perturbative 
calculations using a space of perturbative states does not enjoy the same protection that the corresponding practice 
has in QED. In Section Hill we will establish an isomorphism between the gauge-invariant states that implement the 
non-Abelian Gauss's law and perturbative states. This isomorphism enables us to substitute "standard" calculations 
with perturbative states for prohibitively difficult ones with gauge- invariant states. By this means, we provide for 
QCD a substitution rule, similar to the one available in QED, that permits the use of perturbative Fock states in 
scattering calculations with the assurance that the results of these calculations will agree with what would have been 
obtained if gauge-invariant operators and states had been used. 

Another difference between QCD and QED is related to the fact that states that obey the condition 

g'^(r)|*)=0 (48) 

are not normalizable. We can see this easily by constructing, for example, the commutator of G^(j) and an integral 
operator I = j (ir'A^(r')xj (r') where Xj(r') is an arbitrary c-number- valued function. Since 

[I,g\v)\=^^,x^{v). and {^[l ,g\v)\\^) = ^^,x^{v)m^) , 

and since Q^{y) is hermitian so that (5'|tJ'^(r) = as well as tj®(r)|4') = 0, this leads to (^'1^') = 0, in contradiction 
to the assumption that |^) is normalizable. This argument is a simple extension of one that was applied to the 
Fermi subsidiary condition for QED. ,,2^] In the case of QED, however, this difficulty can be remedied because the 
non-normalizability of the states that are annihilated by the Abelian Gauss's law operator is entirely caused by the 
unobservable longitudinal nonpropagating photon "ghost" modes, which coincide exactly with the pure gauge degrees 
of freedom, and which can be kept separate from the gauge-invariant transversely polarized propagating photons in a 
variety of ways. In QCD, however, transverse modes can be pure gauge, and we do not know of a similarly satisfactory 
resolution of the non-normalizability of the state vectors that satisfy Eq. H48|l . [27l The previously-mentioned 
isomorphism, which will be demonstrated in Section llTTI mitigates this difficulty by establishing an equivalence between 
matrix elements evaluated with gauge-invariant states that are not normalizable, and corresponding ones evaluated 
with perturbative states. 



B. Relation to QCD in the Coulomb gauge. 

Unlike the Weyl-gauge formulation of QCD, in which one can simply set = and impose canonical quantization 
rules on the remaining fields, |2^|33| the quantization of Coulomb-gauge QCD requires that constraints be exphcitly 
taken into account. In constrained quantization — one procedure for implementing consistency with constraints — 
this consistency is maintained by means of the so-called "Dirac-brackets" , which replace the canonical equal-time 
commutation rules. When constrained quantization, such as the Dirac-Bergmann procedure. |3l| is applied to the 
Coulomb gauge, the generator of infinitesimal gauge transformations becomes a constraint; it then must commute 
with all fields, which therefore are invariant to small gauge transformations. Under these circumstances, the gauge 
field would automatically be invariant to small gauge transformations, although it might have discrete numbers of 
gauge copies. 

However, carrying out the constrained quantization of QCD in the Coulomb gauge is problematical; one impediment 
stems from operator-ordering ambiguities of multilinear operator products. For example, in constrained quantization. 
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the matrix of constraint commutators must be inverted. There are noncommuting operators in that matrix, and it is at 
best problematical to keep track of operator order in the process of finding this inverse. As a result, the Dirac brackets 
of some operators are not unambiguously specified. Because of the difficulties associated with the quantization of 
QCD in the Coulomb gauge, a number of workers have avoided the direct quantization of Coulomb-gauge QCD, and 
have proceeded by treating the Aq — gauge fields as a set of Cartesian coordinates and the Coulomb-gauge fields as 
a set of curvilinear coordinates, and have transformed from the former to the latter by using the familiar apparatus 
for such coordinate transformations. 133. 1331133. 1351 



In our work, we transform from the Weyl gauge to a representation in terms of gauge-invariant operator-valued 
fields. Our purpose is to implement gauge invariance, not to carry out a gauge transformation. We do not impose 
transversality on the the gauge- invariant Aq^^; in our work, Aq^^ is transverse, but the transversality is not imposed 
as a condition — it emerges as a consequence of its gauge invariance. And the Gauss's law operator CJ" does not 
vanish identically; in our work. Gauss's law is a condition on a set of states (the implementation of Gauss's law by 
imposing it on a set of states is also discussed in Refs. [13, 13 HE 113 ) ■ 

Because our formulation of QCD in terms of gauge-invariant fields differs significantly from those whose purpose is to 
construct the QCD Hamiltonian in the Coulomb gauge, it is of interest to inquire how closely the resulting Hamiltonians 
resemble each other. In order to examine this question further, we will make some additional transformations of Hq\ 
that assume that the Hamiltonian acts only on states that implement Gauss's law. When Hci appears in a matrix 
element between two states I^I'q.) and (v!//?! that obey G'^i'x.) I^'q) = and (5'^|t/^(x) = 0, further transformations that 
eliminate the longitudinal component of IIq ^ are possible. For the case that Hq^^ appears adjacent to and directly to 
the left of such a state we can make the replacement 



naM^> = -|iJ5(Gi)l^> (49) 



and. therefore, also 



-/o^G,)!*) = ^r'^^H (n^u + n^h) l*> = {^o(g,) - ^r'^^Lf^^iG,)} I*) > (so) 

where Jo(Gi) defined as Jo^Q\^=9f'^^'^J^ai^'Qli ■ Eq- can be iterated, leading to 



•^O(GI) 



n=0 



where « indicates that the replacement is valid only when the operators act on states jvE') that implement Gauss's 
law. When Jq(q) stands directly to the right of {'^\ states, we can similarly make the replacement 

•^o(G,) « - E(-l)"5"r'" (T(?Vo(Gi)) (52) 



where 



= (((«y)(r)) %^ <t)«) <l2)(r) j ^ <!r)(r)j (53) 

and where the arrows indicate that differentiation is applied to the left. Similarly, IIq j(r) and ngp-(r) can be expressed 
as 

n^i, « n^^^. + % ( £(-i)".9".r'"' (%4(G,)) ] (54) 

and 

« n^^/ + % [Y,{^lTg-f'^ (^(n)^o(G,))) (55) 

\n=0 / 

respectively. We can combine Eqs. (|2t)|l with (|51|l and H52|) to obtain 

^o(Gi)(y) « d^j dxI?''«(y,x)Jo"T,)(x) (56) 
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and 



Eqs. H54|) and (|55|l can be expressed as 

n^„(y) « n^u(y) - 9, J dxi?'«^(y,x)Jo"T,)(x) 



and 



(57) 



(58) 



(59) 



respectively, where >/o(go(^) represents the hermitian adjoint of Jp(^|^(x). 

We can define an "effective" Hamiltonian (ffGi)pj^yg, which is obtained by making the replacements described by 

Eqs. H56() - (ISOfl l in Hq\ and excluding Hg, since the latter will not contribute to any matrix elements in the physical 
space in which Gauss's law is implemented. With these replacements, we obtain 

(^Gi)phys = /^r [ in^,V(r)ngT(r) + lF^,,^{r)F^,^^{r) + ^Hr) (/3m - za.a,) V(r)] - j dvj^{v)Al,^{v) 

-\ I drd^dy (jo^(x)) + JoYJ)(x)) V^^ (x, r)9^P-(r, y) (jS(y)) + J^l^iv)) ■ (60) 



(-f^Gi)phys identical to Hq\. But (TJcOpj^yg can substitute for Hq\ as the generator of time-evolution when we 

embed the theory within a space of states that satisfy the non-Abehan Gauss's law, t/"(x)|^i,) = 0. Because 
G"'{x) is hermitian, the same state I'^v) that obeys Eq. (|^ also obeys (^'i/|C/°(x) = 0. Eq. 1^7)1 demonstrates that 
when Hq\ appears in any "allowed" matrix element, IlQ j and Jq(^qi^ always are situated where they abut a "ket" 
state vector \'^a) to their right; and IIq^,- and Jq/gi) ^I'^^iys are situated where they abut a "bra" state vector (^'/jj 
to their left. Since Hq] will always be bracketed between two states and \'^a) that implement Gauss's law, IIq ^ 
and IlQ^j can be replaced by their "soft" equivalents shown in Eqs. H58|l and (|59|l respectively; and Jg^gi) and t/Q^^gi) 
can similarly be replaced as shown in Eqs. (|56() and H57|l respectively. (^^Gi)pJiys '^^^ therefore always be substituted 
for Hq\ in matrix elements, as long as attention is paid to the need to restrict the space of state vectors to those that 



implement Gauss's law. For example, exp{—iHQ\t)\'ii a) can be replaced by exp *(^Gi)phys^ 
be required to project onto states that implement Gauss's law, as shown by 



exp(-i(7?Gi)OI*a) = |*.)(*.|exp(-i(i?Gi)t)|*a), 



I^'q), since both will 



(61) 



and 



{■^^\exp{-i{Ha)t)\^a) = - ii(*^|i/Gl)|*a) + • • • + 

(vI/^J7?gi)|*m.>(*mJ^gi)|*m3 



(^'.li^Gl)!^' 



(*M„-J^Gl)|*o) 



(62) 



Each matrix element ('J'^Ji^GOI^Mj) '^^ ^Iq. (|^ can be replaced by (vl/^. |(iJGi)p}jyg|*^j), so that cxp(— i(7?Gi)^)|*a) 
can safely be replaced by exp (^^*(^Gi)p}iygi^ I^q)- The time evolution imposed by Hq\ on a state vector |5'q) for 
which ^^^(x)!^'^) = takes place entirely within the space of states that implement Gauss's law. In the case of a state 
vector \x) for which 5'^(x)|x) = \x') where \x') is nonvanishing, 



(Xl g^(x) ,exp(-ii?Git) |*a) = {x'\exp{-tHama) - 



(63) 



because ^"^(x) and Hq\ commute. This requires the part of x that fails to implement Gauss's law to be orthogonal 
to exp(— i_ffGii)|^Q:)- The only limitation on the validity of this argument is the non-normalizability of the states 
that implement Gauss's law, which complicates the algebraic properties of the {l^a)} vector space. Nevertheless, 
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Eqs. (|61() - (|63|l show that we can restrict the space in which time evolution takes place to state vectors that implement 
Gauss's law without compromising the unitarity of the time-evolved \'i'a(t)) or of the S'-matrix evaluated with such 
states. These considerations are also instrumental in allowing us to replace exp{—iHQ\t) with exp(— i(ifGi)phys^)- 
and (iJcOpiiys t'oth commute with t/"(x) for all values of a, so that 

^'^(x)exp(-z(i7Gi)phys^)l*o) = cxp(-z(7?Gi)phys<)^"(x)|«'o) = (64) 

as well as 

g'^(x) exp(-i(i?Gi)t)|*a> - exp(-z(HGi)i)a''(x)|*„> = . (65) 

The state vectors eyip{—i{HQ\)t)\'i/ a) and exp(— i(i/Gi)pJiys^)l'^«) therefore are gauge-invariant and implement Gauss's 
law just as \'^a) does. 

In comparing (i/Gi)p}^yg with expressions for the Coulomb-gauge Hamiltonian in the literature, we note that the 
only significant difference between (-ffGi)phys ^he Coulomb-gauge Hamiltonian reported in Ref. (s^ is that n^,'^.^, 
the hermitian adjoint of the transverse gauge-invariant chromoelectric field, appears in Eq. H60() where the expression 
J-^n^Lj appears in Ref. where J = det[a,-A]. We will prove in Appendix IbI that 



KJ,^ = J-'^lJjJ , (66) 



by using Eq. (|ll(l and the identity 



(67) 



where the trace in Eq. H67|l extends to the coordinates and the color indices. With this demonstration, we see that 
Eq. (|60|l and the Coulomb-gauge Hamiltonian described in Eq. (4.65) in Ref. [s^l are identical. It is also of interest to 
compare Eq. (|60(l with the Coulomb-gauge Hamiltonian in Ref. 11] as well as in the work of a number of other authors 
who used the same form of the Hamiltonian. The Hamiltonian in Ref. differs from the Hamiltonian described by 
Eq. (4.65) in Ref. [3^ in the fact that H^^- rather than H^/"/ appears in Ref. [13 in place of J'^UlJ^J in Ref. 
there is also the trivial difference that Ref. |^ deals with "pure glue" QCD so that the quark field is not included. 

This discrepancy raises the question of the hermiticity of the operator-valued transverse gauge-invariant chromo- 
electric field ^qJj, which is of considerable importance for determining the dynamical effects of (_ffGi)pi^ys- One way 
of addressing this question is to use Eq. Hll() and Eq. (65) in Ref. to obtain 

K^Jiy)~'^GAy)^[^]iy),Rb,iy)]=^9f''''^V^'^{y,y) (68) 

where the partial derivative acts on only the first y argument in I?^''(y,y). We might have expected that the 

transverse parts of HQ^ (y) and HQ^ (y) would be identical since any functionals of the form {Sij ^)dj£,{y) would 

necessarily vanish. Such a conclusion would not, however, be correct in this case, because in -^'D'^^{y, y), the partial 

derivative differentiates only the first y in I?'''*(y,y). We can make use of Eq. 1)38(1 and the fact that f^'^'^Shc = to 
express Eq. as 

n^L(y) - ^Uy) - ^9^t''f-j d^^^ (i^) 4,.(z) Ap..(,,y) (69) 

and we can extract the transverse parts to obtain 



Eq. iffnji makes it clear that Hg|'^^(y) — Hg|'^ (y) is not the transverse projection of a gradient and therefore cannot be 
presumed to vanish. 
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Equally compelling evidence that ^qIj is not identical to its hermitian adjoint is provided by the observation that 
ngi"^^(x) ,ng/"/(y)j and [ngT(x) , (y)j differ. The latter vanishes, as is shown by Eqs. (gOl-lllII. 



the commutators 

However, use of Eq. Ijlll) and the commutation rules for the underlying Weyl-gauge fields lead to 



2 rhca 



9 f 



9-^ 

92 



92 



d 
dyi 



I?'"^(y,x)gf^I?^P(x,y) 



(71) 



and an alternate derivation based on Eqs. H42|) and (|68|l confirms that result. Similarly to what we observed in 
connection with Eq. H68() . the derivatives and g|- each differentiate part, but not all of the y and x dependence, 

respectively, of the product I?'^'*(y, x)I?'^p(x, y) in Eq. (|71|l . The transverse projections of ^I?'"'(y, x) g|^I?'^P(x, y) 



therefore will not vanish, and 



n^uW,n^;'/(y)J ^ 0. Since [n^T(x),n^T(y) 

n^i'^ (y) and nQ'^^(y) cannot be identical. 



and 



ns,"^.(x),n^^t(y) 



differ 



III. ISOMORPHISM AND ITS IMPLICATION FOR THE SCATTERING AMPLITUDE. 

In the preceding sections we have obtained a description of QCD that took the Weyl-gauge formulation as its point 
of departure, and arrived at a Hamiltonian in which all operator-valued fields — the gauge field, the chromoclcctric 
field, as well as the quark field — are gauge- invariant, and only the transverse components of the chromoelectric fields 
appear in the Hamiltonian (^^Gi)phys- necessary, however, to restrict use of this Hamiltonian to a space in 

which all state vectors implement the non-Abelian Gauss's law; and these state vectors are complicated constructions 
that are not easy to use. In this section, we will show how isomorphisms can be established that enable us to identify 
(^Gi)p}iyg with a Hamiltonian that can be used in a space of ordinary, conventional perturbative states. 

To review the relation between gauge-invariant and perturbative states: In Ref. 0, a set of states was constructed 
in the form 

1^.) = ^l*^.) (72) 

where the operator-valued 'J was given as 
with 

= / dr,---dr^:^^in)---^^ir^)nf^]{r,)--^^ (74) 

designates one of a set of states that is annihilated by djU^. These \4>i) states — the so-called "Fermi" states — 
are related to "standard" perturbative states \pi) by 

\^^) = S|p,) ; (75) 

S was given in Ref. [sl], where it was also shown that 9jHj(r) S|pi) = 0, where |pi) designates one of a set of "standard" 
perturbative states annihilated by all annihilation operators for fermion and transverse gauge field excitations. This 
set of perturbative states will be described more fully later in this section, and will turn out to be identical to 
perturbative states in QED, except for the fact that the gluon operators carry a Lie group index, while the photons do 
not. Since 9jHj annihilates any {ipi) state, we can see that, in states, the negative chromoelectric field Tl'^^^^iji) 

in ^I* can be replaced by its transverse part H^^^^^j^(r£), because the longitudinal parts vanish when acting on a 
state. Furthermore, in Eq. (|74|) . every transverse H^'^j^(r^) is integrated with an ^^^^^^^(r^) in each variable, r^, and 
only the transverse components A'j}^^-^'^ (rg) will survive this integration in the j^'^) states, which become 

■n n 

1*^) = Y.^ dr,---dr^ AlWY ir,) ■ ■ ■ Af^f ir,,)uf,Y (r,) ■ ■ ■ Hg:)/(r„)|0.) . (76) 

n=l ■ 
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In Ref. fl^, it was shown that 

Ag,/r)* =vI/Af(r)|0,). (77) 
In Appendix IdI we will use Eq. ()76|l to show that 

n^T(r)vI/|0,) = vI/nf (r)|0,). (78) 

Since the Hamiltonian (iJcOpjjyg consists of transverse fields only, Eqs. H77|l and H78|l afford us an opportunity to shift 
(^^Gi)phys from the left-hand side of 5* to the right, with a concomitant substitution of transverse Weyl-gauge fields 
for the corresponding gauge-invariant fields. The one impediment to this process is that n^i*^^, the hermitian adjoint 
ofYlolj, also appears in (^^cOphys' ^^"^ o^^y applies to fl^/^ and not to H^i'^^ To remove that impediment, 

we use Eq. (jnni) to substitute J-^n^(;JjJ for Hg/"/, and express (^Gi)p}iys as 

(HGi)phys = / rfr {i^-in^T (r):7ngT(r) + ii^G^.W^^GV,, W + V^^(r) ipm~taA)Hr) - 
jT(r)A^„(r)} - i I drdxdy (jo^(x) + JoYJ)(x)) T^^ (x,r)a2p-(r,y) (j-(y) + Jo^(T,)(y)) 

(79) 

with 

'^o(Go ~ gP^'^J ^ng/^j'AQj . 

We can define a "hermitized" transverse gauge-invariant negative chromoelectric field V^~^ 

V'^\v)^J-'^nlJ^{v)j'^. (80) 
As can be seen from Eq. (|66|l . V^~^ is hermitian, since 

V]^\v) ^ J-m'^,^{v)J-'^ = J-'^ {jTi'^,^{v)J-') = V',\v) . (81) 

An important consideration for this argument is the fact that is hermitian, which is proven in Appendix El In 
the same appendix, we also prove that the canonical commutation relations between H^^'s and ^q^ 's and that among 
n^^'s remain unmodified with H^^'s replaced by VqI^'s. We then find that 

n^u(r) = jiv^^{r)J-i and Il'J^^r) = J-iv^^{r)ji . (82) 

Eq. H82(l transforms from the nonhermitian fl^i'^ and to the hermitian Vj~^ (not, however, unitarily, since J^^ is 

hermitian and not the hermitian adjoint of J^^^). Transformations of this kind have previously been used by other 
workers; [s^. |3^ It would be possible to make a compensating transformation on the states, but we prefer to leave 
the states untransformed and to extract 

[H]o = jdv[ ^r^^{r)V^^{r) + iFc^H, Wi^GH, W + ^^(r) (/3m - laA) V'(r)] , (83) 

from Eq. (|79() in order to obtain a non-interacting part of (ifcOpi^ys that consists of hermitian gauge- invariant fields 
and that can define interaction picture operators. As we will show in Appendix IeI this process leads to the expression 

(^Gl)phys = [H]0 + [H]l + [H]2 (84) 

where 

[H]i = J dr {gr^-a,Ag,^.(r)A^,,(r)A^,^.(r) + igV"''V"'''^'^^,.(r)A^„(r)4,,(r)Aa,,(r) 

- jT(r)Ag„(r)} - i J dvd^dy (jo^(x) + Jg[G,)(x)) V^^ {^,r)d'V^^{r,y) (j-(y) + JS[G,)(y)) (85) 
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and 



U + V + ^J dvd^dy |»kg(x) V (x, r)d'V -(r, y) (j-(y) + Jj^lci) iv 



)) 



j;;(x)+j;;^C,)(x)j P"'' (x,r)9^I?-(r,y)ikEj(y) + k[5(x)P'«^ (x, r)a^I?-(r, y)kEj(y) 

where and V as well as kQ(x) and Jq (g,) are defined in Appendix lElin Eqs. HE7|) . IjElip . (|E16p . and ljE15() respectively. 
[H]o, [H]i, and [H]2 are hermitian, and all consist entirely of gauge-invariant, hermitian, transverse gauge fields and 
gauge invariant quark fields, which all obey "standard" commutation rules. Since 7^j'^(y) and ^Q]j{y) have the same 
commutator with (x) , Eq. (|42|l also determines the commutation rule 

[^Gi,(y),^GHW] = [V^^{y),V^^{^)] =0, 



Q2 



(87) 



The sum [H]o + [H]i is identical in form to the Coulomb-gauge QCD Hamiltonian used by Gribov 0,0], as well 
as by numerous other authors who have followed him in using this Hamiltonian. [H]2 consists of additional terms 
that are required because the transverse gauge-invariant negative chromoelectric field ^qJj is not hermitian. The 

elimination of and ^qIj^ in favor of the hermitian Vj^ is essential for the establishment of the isomorphism 
between (iJcOpi^yg '^^'^ ^ Hamiltonian that can be used in a Fock space of perturbative states. We now proceed to 
the demonstration of this isomorphism. 

Since both j(r) and 7'j'''^(r) are hermitian and obey the commutation rule displayed in Eq. H87() . we can represent 
them as 



A^, ,(r) = J2 ^ K(k)e'^''-- + <Hk)e-^^-] 



(88) 



and 



k, n 

where n is summed over two transverse helicity modes and 

<(k),a,^^(q)] ^Sru6a,bSi.,^ and [<(k) , a,^(q)] = [a^^k] ,a';\q} 
Eqs. (|88|l and H89I) can be inverted, leading to 



<(k) 



e."(k)/dr( A^,,,(r)-f-7'r'(r))e 



-zkr 



and 



(89) 



(90) 



(91) 



(92) 



Eqs. (|91|l and H92|l show that afj(k) and a^t(k) are gauge-invariant and commute with the Gauss's law operator tj"(r). 
Eqs. (|77|l . (|78|l and H8()(l demonstrate that any functional T (^gu '^j^) ^^^^ have the transformation property 



leading to 



(r) e 



-ikr 



\4>i). 



(93) 



(94) 
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and 

^ ' ' ^ ^^^{r)]e-^^-}\cl^,), (95) 



so that the isomorphism estabhshed in Eq. H93II between the gauge- invariant fields j, Vj'^ and the gauge-dependent 
Weyl-gauge fields A"^, respectively is transferred to a similar relation between the gauge-invariant creation and 
annihilation operators for transverse gluons, a^^(k) and afj(k), and corresponding "standard" perturbative creation 
and annihilation operators a^^^ (k) and (k) . We can proceed by using the standard representation for the transverse 
part of the Weyl-gauge fields, 



At^v) = K(k)e^'^-"- + <,t(k)e-''--] (96) 

' — V 2k 

k, n 



and 



nfir) = -^Y.^^'\^)^[\ K(k)e'''-- - a'i,\\.)e-^'^-] , (97) 

k, n 



which demonstrate that 



a'J{-k)J—^^\cl),) = J—^^a'J{-k)\cf>,) and 

^(k)^-^*!^.) = J-Ha'f,{kM,) . (98) 

Any a^(k) will annihilate the gauge-invariant vacuum state j/^^5'S|0), because the transverse excitation operators 
afj(k) and a^^(k) trivially commute with S. 

At this point, we can establish an isomorphism between two Fock spaces: The "standard" Weyl-gauge Fock space 
consists of 



(k)|0) (99) 



|ki . . . k, . . . k^r) = K K\t(k,) . . . a'^;{k,) ■ ■ ■ aZH^N)] |0) ; (100) 

with K the normalization constant and the gauge- invariant states that implement the non-Abclian Gauss's law can 
be represented as 

\k) = ljalHm~^mO) (101) 



|ki . . . k, . . . kAT) = f [<t(ki) . . . <;t(kj . . . aZH^N)] J-imO) (102) 

where |0) designates the perturbative vacuum annihilated by a^(k) as well as by the annihilation operators for 
quarks and antiquarks, ^p^s and qp^s respectively. The additional normalization constant must be introduced to 
compensate for the fact that |Cp = | J'~^^'S|0)p = (0|S*\I'*J7'~^^'S|0), which formally is a universal positive constant, 
is not finite; and the state JT'"^ VPSjO) is not normalizable. However, once C is introduced, the |ki • • • ki • • • kA?) states 
form a satisfactory Fock space that is gauge-invariant as well as isomorphic to the space of |ki • • • k^ • • • k^r) states. 
We can now use Eqs. H88|l . and H89() to express [H]o as 

[H]o = Y. (k) + E Wp,s%,s + <ils%..s) (103) 

k,c p,s 

with the subscript s labeling the color, flavor and herlicity of the quarks. In this form, [H]o can be seen to describe 
the energy of non-interacting gauge-invariant transverse gluons of energy k and quarks and anti-quarks respectively 
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of energy £p = A/m2 + |p|2. We can also define another Hamiltonian, Ti = Tig + Tii + 7^2, in which each component 
part is identical in form to [H]o + [H]i + [H]2 respectively, but with the substitutions 

P^^T^nJT and Al,,-^At^ 

everywhere — including the replacement of Aq^ ^ hy A°;~^ in the inverse Faddeev-Popov operator D °''(x, y) — so that 
Ti is characteristic of the Coulomb gauge, but nevertheless is a functional of transverse Weyl-gauge unconstrained 
fields. For example, Tio is 



Ho^ dr 



inf T(r)nf T(r) + i^;;;.(r)i^'5.(r) + ^p^r)if3m-iaA)ib{r)\ (104) 
where Fj" — djA'^ — diA"-. Using Eqs. and we can express Tio in the form 

Ho = E (l^X (1^) + E ('^P.^'^P.^ + ^l,s^T>,s) ■ (105) 

k,c p,s 

We can then use Eq. to establish that 

[H]o^-5$ = SHo|n) (106) 

as well as 

[H]iJ-5*S|7i) = J-^*SHi|n) (107) 

and 

[H]2^-5$ S|7i) = SH2|n) (108) 

The state vector \n) represents one of the |ki • • • • • • kjy), the "standard" perturbative eigenstates of TYq- 

We can use the relations between Weyl-gauge and gauge-invariant states we have established in the preceding 
discussion to extend the isomorphism we have demonstrated to include scattering transition amplitudes. For this 
purpose, we define 

T^int - Til + H2 and [H]i„t - [H]i + [H]2 . (109) 
The transition amplitude between gauge-invariant states is given by 



T^,, = ±(/|S*Vl,*j7-t } + [H]int7 — Mint ) J'^-^m , (110) 

^ ' (i?.-(i/G,)phys+*^l 



where \i) and |/) each designate one of the \n) states; \i) represents an incident and |/) a final state in a scattering 
process. With the results of the preceding discussion, we can express this as 



= ^{Q\S^^'*J-Hm{f\ (Wint +Hi„t-— -— ^— — -Hint j K) (111) 

y (El - Ha- Hint + «e) J 

where we sum over the complete set of perturbative states The second line of Eq. 1111|) follows from 

i^(/|S**V-^*S|n) = ij(0|aKk/)S**V-^*Sat (k„)|0) = ij(0|S*vI/'^^-ia/(k/)at (kn):^-**S|0) = 

5f^J{\^f - k„)ij(0|S**V-^*S|0) - ij(0|S*vI/V-*at (k„)a/(k/)^-ivl,S|0) (112) 

and the observation that the last term on the second line of Eq. pi2|) vanishes trivially. With the isomorphism of the 
states |ki • • • ki • • • k^v) and |ki • • • ki • • • Icn) that we have established, 

T/,. = T/,, (113) 
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where 

Tf,^ = (/I (n^nt + Hint .^ ^ ^ . . Hint) N) (114) 

1^ (iii - rto - /tint + «ej J 

is a transition amplitude that can be evaluated with Feynman graphs and rules, because it is based on "standard" 
perturbative states that are not required to implement Gauss's law and need not be gauge-invariant. 

In the remainder of this section, we will discuss the relation of our formulation of the scattering transition amplitude 
to approaches to this problem in Coulomb-gauge formulations of QCD. As was pointed out in section IIB, the effective 
Hamiltonian (i?Gi)phys described in Eq. (|60|l is identical to one obtained by Christ and Lee, js^ who treated gauge 
fields as coordinates and applied the apparatus of transformations from Cartesian to curvilinear coordinates to the 
problem of formulating Coulomb-gauge QCD. Here, we will show that Hq\ — the precursor of (J?Gi)phys , described in 
Eq. H2(J|I — is identical in form to the Hamiltonian given in Eq. (6.15) in Ref. [s^, which leads to the Coulomb-gauge 
perturbative rules formulated by Christ and Lee. For this purpose, Hq\ will be expressed in terms of Vj^ and Aq j, 
and then Weyl-ordered. The equivalence of Christ and Lee's results with Schwinger's [s^ was already confirmed in 
Ref. '3^. 

Eq. (|82|) demonstrates that the functional dependence of Hq\ on T^j"^ and is the same as the functional 

dependence of 

H = J^-HaJ--2 (115) 

on HJ^ and Aq,;. H was used by Christ and Lee to generate the path integral representation of the Coulomb 
gauge, and they showed that 

iJ = 7?^+Vi+V2 (116) 

where the superscript W designates Weyl-ordering with respect to Hgi*^ and Aq^. The additional terms Vi and V2 
are given by 

Vi = \9^f'"f''^ J drd^V-'ir, v)d,V-\v, r) (117) 

and 

V2 = iff'/'""/'"" j rfxdydz (5kn5,,5{y - x) + i?jl?'="(y, x) a. 

V--{^, ■L)d''V-\'L, y)) - y) + A2?'"(x, y) 9^) (118) 

where the partial derivative dj to the left of I?°''(r,r) acts only on its first argument. When a partial derivative with 
a left arrow on top appears to the right of I?"^ with two identical arguments, it acts only on its second argument. 

The case of of two identical arguments of T) is understood as 9jl?'^''(x, x) = limy^x gf-2'°''(x, y) and I?"''(x, x) d j = 
limy^x gf~^''''(yj x), where the limit is taken after the partial derivative has been evaluated. This convention will be 
followed consistently in the following discussions. Since the commutator of and Aq ^ is identical to that of Hgi*^ 
and , an equation parallel to Eq. (|116|l . 

i^Gi = + Vi + V2 (119) 

will be proven below. The superscript W designates Weyl-ordering, but in this case with respect to and ^ . 

The parallel structure refers to the fact that, as was pointed out above, Hqx has the same functional dependence on 
and Agij as ¥l has on ^(^j and ^Qj. Since the fermion variables commute with V^'^ and ^Qj, we may drop 
them for the proof of Eq. 1)119(1 ; we will also drop Hg , since it makes no contributions in the space of gauge- invariant 
states. It follows from Eq. ((58() that for a physical state j^*), 

n^,^.(r)|vl/)=-4(r)|vI/) (120) 

with 

4(r) = + y dx9,I?^9(r,x)A(x)n'^(x)]j-^. (121) 
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In terms of the Weyl-ordered chromoelectric field operator of Schwinger, |^ 
we have 

£!>{r) = E'^{r) + A';{r) and 
f}(r)t=£;^(r)-A5(r). 

The Hamiltonian Hq\^ in the absence of the fermion field and without Hg, can be written as 

1 



(122) 



where the kinetic energy 



with 



To evaluate A^(r), we observe that 



dv[E]{v)E]{v)+v{r)] 
«(r) = ~A5(r)A?(r)+[K^(r),A5(r) 



(123) 



(124) 



(125) 



(126) 



£b{r) = ~V'^^(v) + \ [V';^ {v) MiJ)] + j d^dp"i{v,^)D,{^)vf {^)-\ J dxa,I?^«(r,x)A [^''^(x) , ln( J)] (127) 
and 



so that 



and, therefore, that 



vf{^)M{j) 

(128) 



1 ( cb 
2 



£:}(r)+f^^(r)t)=£;)(r) 



With Eqs. (|nn|l and (gHl), this leads to 



In Appendix^ we shall prove that 



\{£]{v)^£]{v)^)^A]{r). 



i j dvv{v) - Vi. 



(129) 
(130) 

(131) 
(132) 



In the form given in Eq. (|124|l with K as described in Eq. p25|) . the effective Hamiltonian {Hq\) is identified with 
that of Schwinger. i32j The next step towards the proof of (|119|l follows from the operator identity given in Ref. [l^ 



- / dvE'^Av)E){r) - - / dr £;,^(r)ii;^(r) 



w 1 

— - I dxdydz 



(133) 

Using the commutation relation (jElOp . we can show that the second term on the right hand side of Eq. (|133|l is the 
same as V2 (the same proof is also given in Ref. [s^) and Eq. pi9|) is established. 
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IV. DISCUSSION 

In this work, we have used earher results 0, 0, @| to express the Weyl-gauge Hamiltonian entirely in terms of 
operator- valued fields that are gauge- invariant as well as path-independent. These gauge- invariant fields have many 
features in common with Coulomb-gauge fields: Their commutation rules agree with those given by Schwinger in his 
Coulomb-gauge formulation of QCD, [3113 except for differences in operator order; these differences can be ascribed 
to the fact that Schwinger imposed Weyl order in his work while we do not make any ad hoc changes in operator order. 
The gauge-invariant gauge field is transverse and hermitian; but the gauge-invariant chromoelectric field is neither 
transverse nor hermitian. Even the transverse part of the gauge-invariant chromoelectric field is not hermitian. That 
fact is important for relating the Hamiltonian we obtained in Eq. H84() with those given by Gribov, Tl'l Schwinger, l3a| 
and Christ and Lee. 

The relation between the Coulomb-gauge Hamiltonian for QCD and the Weyl-gauge Hamiltonian expressed in 
terms of gauge-invariant fields closely parallels the relation between the two corresponding QED Hamiltonians. The 
Weyl-gauge Hamiltonian for QCD is represented entirely in terms of gauge-invariant fields in Eqs. (|20|l and H24|l . 
When formulated in terms of gauge-invariant fields, QCD must be embedded in a space of gauge-invariant states 
that obey the non-Abelian Gauss's law. Within such a space of gauge-invariant states, further transformation of the 
QCD Hamiltonian we have constructed can be effected. Thus transformed, the Hamiltonian consists of two parts. 
One part, (-ffcOpJiyg — displayed in Eq. H6()|l — is identical to the Coulomb-gauge Hamiltonian. It is a functional of 
transverse gauge- invariant chromoelectric fields, gauge-invariant gauge fields (which are inherently transverse), as well 
as gauge-invariant quark fields. The other part, Hg — displayed in Eq. (|24ll — makes only vanishing contributions to 
matrix elements within the space of gauge-invariant states that is required for the Hamiltonian to act consistently as 
the time-evolution operator. Hg does affect the field equations and "remembers" that the formulation is for the Weyl, 
and not the Coulomb gauge. This situation is precisely the same as in QED, in which the Weyl-gauge Hamiltonian, 
expressed in terms of the gauge-invariant field (in that case, simply the transverse part of the gauge field), is the sum 
of two terms, given in Eqs. H45|l and H46|l: the former is the Coulomb-gauge Hamiltonian, and the latter makes only 
vanishing contributions to matrix elements within the space of gauge-invariant states, but is necessary for reproducing 
the Euler-Lagrange equations for Weyl-gauge QED. 

In spite of the similarity between QCD and QED in the relation between the Weyl and Coulomb gauges summarized 
in the preceding paragraph, there is an important difference between the gauge-invariant states for the two theories: 
Gauge- invariant and perturbative states in QED are unitarily equivalent; and in a Hamiltonian formulation, this 
unitary equivalence permits us to use perturbative states in evaluating scatteri ng amp litudes in QED in algebraic 
and covariant gauges without compromising the implementation of Gauss's law. |22ll23l| But there can be no unitary 
equivalence between gauge-invariant states and perturbative states in QCD. And the gauge-invariant states in QCD 
are complicated, not normalizable, and very cumbersome to use. In order to make effective use of the Weyl-gauge 
QCD Hamiltonian represented in terms of gauge-invariant fields, some relation is required that allows us to circumvent 
the absence of the unitary equivalence between gauge-invariant and perturbative states that afflicts non-Abelian gauge 
theories. In Section Hill we establish such a relation in the form of an isomorphism that enables us to consistently 
carry out calculations in QCD with an equivalent Hamiltonian that is a functional of the original gauge-dependent 
Weyl-gauge fields and that is used with standard perturbative states. In the case of QCD, this isomorphism has 
been demonstrated for the Weyl gauge only. An extension to a somewhat larger class of algebraic gauges defined 
by Aq 7A3 = with 7>0 should not be difficult; ^3 but, in contrast to QED, there is no indication that further 
extensions — to covariant gauges, for example — are possible. Finally, in Section IIIII we show that the effective 
Hamiltonian (/fcOphys — and therefore also H — Ho + Hi + H2 — can be expressed in appropriately Weyl-ordered 
forms and shown to be equivalent to results obtained by Schwinger js^ and by Christ and Lee |32l | . The Hamiltonian 
used by Gribov in Ref. [ij is equivalent to only H = Ho + Hi- H2 does not appear in that work, because the 
nonhermiticity of the transverse chromoelectric field was not taken into account. 
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APPENDIX A 

In this section, we will prove Eqs. (|31|l and H3U|) . We use Eqs. (|34|l - (|36f) and expand the product 'D°-^{y, x) 9--D''''(x) 
in Eq. H31|) as a series in powers of g, and observe that 0{n) terms originate from I?"^'^(y,x) (^hb (x)^ a-nd from 

(y, x) ^—.9/'"''' di • For example, the first part of the ri = term of Eq. H31() originates from the 6hb 

part of d-D^^ (x) and is 



, , , , 5hb (X) = a^(x) = -5„fc5(x - y) (Al) 

47r|y-x| V 7 M7r|x-y| 

and the second part of the n = term of Eq. H31|) . which stems from the —5/'"''' di Aq ,. in d-D^'^ (x), is 

j;;,^ (-.,/'•" a <„(x)) = -i,r"ft,.,3;5i3^A.,.(x) . (A2) 

The first part of the n — \ term of Eq. H31|) 



(x) (A3) 



exactly cancels the second part of the n = term, and this pattern of cancellation can easily be seen to hold in general 
— the first part of the n + 1 term cancelling the second part of the n-th term. For the general term, 

vili) (y, x)5,, a^x)= . . v/ ^^^^^^^Gi^, (^(1)) X 



47r|z(i) -z(2)| '^i' ' 47r|z(n-i) -x| 

and {I?f„^i)(y,x)}{-5/''9^ a. A^h(x)} is 



47r|z(i) - z(2)| ' "9z(2),2 47r|z(n-i) - x| 

so that, relabeling dummy indices /i— >s„_i and q — > Q!„, we obtain 

2?^„^)(y,x) 9V) -2?^,ti)(y,x).g/''^'' 9, ^^h(x) = 

and a consistent pattern of cancellations is established with 5ab5{y — x) remaining as the only surviving term in 
pa/i^-y^^-j Q.jjhb (^^y similar argument can be used to demonstrate Eq. (|3(J|I . 

APPENDIX B 

In this section, we shall prove the identity given in Eq. (|66|l . By definition — Eq. (|11|) — we have 

^-in^,^.(x)^ = n^,^.(x) + /j'(x) (Bi) 

with 

= j-^ [n^i,(x) = RbcW'^ [n,^(x) (B2) 
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SO that 

= -ii?6,(x) Idyl dzl?™"(y,z)^^^9-l?"'"(z)5(z-y) 
= *.g/'""i?ftc(x) I dy I dzl?'""(y,z)^^ii^^^j(z-y) 

= -.9/'"™i?tc(x) I dy I dzl?™"(y, z) [n^^(x), 4,,(z)] ^8{^ - y) (B3) 
where we have used Eq. (|^ . Substituting Eq. and using 

n^^(x)==i?,,(x)n^i^-(x), (B4) 

we obtain that 

= *5/''^'^^P"^(x,x) + (B5) 

where the gradient acts only on the first argument of T)"-'^ and the longitudinal term comes from the second term of 
Eq. 

/(x) = I dy Ap-(y, y)I?™»(x, y) i^" (y) (B6) 

with acting on the first argument of I?''^(y, y). Comparing the transverse part of Eq. ljB5|l with that of Eq. JHEJ, 
Eq. Hfi6|l is proved. 

APPENDIX C 

To prove the hermiticity of the Faddeev-Popov determinant J as an operator in the Hilbert space of states, we 
recall the criterion that an operator is hermitian if its expectation values with respect to all states are real. In the 
coordinate representation of states for which Agi is diagonalized and corresponds to a c-number field configuration, 
the expectation value of an operator which is a functional of the operator Agi is equal to the same functional of the 
c-number field configuration Agi. For each c-number field configuration, the Faddeev-Popov operator, djDj, with Dj 
denoting the covariant derivative, Dj^ = 5°-^dj — gf^^'^A^^^, becomes an operator with respect to space coordinates 
and group indices. We have 



and 



d] = -d, (CI) 



b] = -D„ (C2) 

with the dagger referring to space coordinates and group indices. Therefore 

id,D,)^ ^D,d,^d,D„ (C3) 

where the last step follows from the transversality of Aq\. Therefore the Faddeev-Popov operator is hermitian with 
space coordinates and group indices for any field configuration. Its determinant, must be real. The hermiticity 
of in the Hilbert space of states is established according to our criterion, and the hermiticity oi J^^ is an obvious 
corollary. 

To derive the commutation relations among 7'j(x)'s and j4q^ (x)'s, we notice that 

Pj(x) = n^T.(x) + 1 [nJ,T (x). In J] (C4) 
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with the second term a functional of (x) only. Then we have 



Furthermore 



where the Jacobian identity 



+i [ngT(x), [u'S,{y)Mj]] + \ [[ngT(x),in:7] ,n^"i;.(y)] 
= [ngT(x),n^T(y)] 



[ngT(x), [n^T (y),in^]] + [[ngT(x),in:7] ,n^"];(y)] 
= -[in^, [ngT(x),n^T(y)]] 

has been employed. Therefore the commutation relations among ^^ (xj's and Aqj (x)'s remain canonical. 



APPENDIX D 

In this section, we shall prove Eq. (|78|) . Using Eq. (|76|) . we define 



from which it follows that 



This leads to 



< = j dv^...dv^A'^;'{v,)...AYJ{v^)ii'^;'{v,)...iiYj{v,,), 



n^(x)*T = vI/Tn«(x)+z / dy n^(x),^f(y) M'Tn5(y 



The commutator involved can be calculated from the relation 

lb 



which implies that 



n^(x),^ T(y) 



(C5) 



(C6) 



(C7) 



The functional derivative was calculated in Ref. |32j| and the commutator in Ref. 5], and the result can 
deduced from Eq. (|^ with the aid of Eq. ((TT|l , which gives rise to 

[n^x) , ^^,^.(y)] = ^ (^RM6,,S{k - y) + i?,,(x) Ap/^(x, y) i5f (y) 

Substituting this into Eq. ljD3|) and using Eq. (|11|) . we find that 

ng,(x)*|0 (x)|</> > dy2?'^'(x,y) of (y)*nf (y)|0 > 

Taking the transverse part of both side, we end up with 

ngT(x)vi/|^>=vi/nf(x)|0>. 

The identity Eq. ((THJ is proved. 



(Dl) 

(D2) 

(D3) 
(D4) 

(D5) 
also be 

(D6) 

(D7) 
(D8) 
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APPENDIX E 



In this Appendix, we will show how to obtain Eq. (|84|l from Eq. (|79|l . In order to obtain the bilinear product 
'Pj~^{r)Vj~^{r) for inclusion in a non-interacting part of (-ffGi)phys that can define interaction picture operators, we 
now express J^iv!!'^ {r)J'P!l'^ {r)J-i 



J-^V^ ' ir)JV!; ' ir)J—2 = V] ' ir)r'> ' (r) + V!; ' ir)J^ ' (r) , ^ 



r'^{r),J--2 J^r'^ir)- V'^{v),J-^ J V'^{r),J-^ 



and make use of the formula 



(El) 



(E2) 



We also observe that the commutator of V^~^{y) and any functional of ylQj(r') commutes with any other functional 
of AQj(r'), and that, in fact. 



ab I ^ij 



Q2 



'^(x-y) 



(E3) 



With these observations, we obtain 

J-ir!l^ir)JV^^{v)J-i = Vy{v)V'^^{v) - \ {Vf(v) , [^'^^(r) , ln(^)] ] 

Eqs. ^ and show that 

With Eq. (|38|l . this can be rewritten in the form 

[n^^,(y) ,ln(^)] = ^9'r''f''^M J A (^^^^) A^,,(z) Ap^''(z,y) 

g(y)a(y) 



(E4) 
(E5) 

(E6) 

where 6j^.{y) = [^Sjk — g"*" j S{y). In this form, it is clear that, to leading order, IIg|^j-(y) ,ln(J') is a term 

and that the limit x ^ y has already been carried out. We can use Eqs. (|E4() and (|E6|) to obtain an expression for 
U = —jfdr [Vj~''{r) , In(J')] [7^j'^(r) , In(J')] , which becomes an interaction term in (iJcOpj^yg, given by 

d 



U 



1^ g4: j^bdh j5ds jhd h' j5' d s 



Jdydzdz' (^Sl„{y)Sl„,{yy 



Av'''{z,y)— , , 
dzi dyk' V47r|y 



dyk V47r|y - z| 



Similarly, from Eq. IjESfl . we see that 



in which we represent 7^j'^(y) as limr^y 7^j'^(r) so that 

[^r(y) ' [^'^(y) ,ln(^)] ] = Imi 5j,(y)^liin A [p^^r] ,V'\y,^) 

Using Eq. (gH), we obtain 

[V^^(r) ,V^^{y,^)^ -gf"' JdzV^%y,z) [T^I^M , Ag, ,(z)] 2?*''(z, x) 



(E7) 



(E8) 



(E9) 



(5(z-r)9^^^X'*''(z,x) 



(ElO) 
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After integration over all of space, J dy [Pj'^(y) , [^^^^(y) , In(J')] ] becomes another interaction term in (^^'cOphys' 
given by 



(Ell) 



V is singular since the leading terms in I?'**(y,r) and T>^^{v,y), {—6ds (47r|y — r|)~^ and —6th (47r|y — r|)~^ respec- 
tively), are not eliminated by the structure constants in V. Christ and Lee called attention to such singularities in 
their work, and conjectured that they might be useful in cancelling unresolved divergences in Coulomb-gauge 
QCD. The same remark applies to V. We continue by eliminating the nonhermitian chromoelectric fields from JI^^q^-^ 

and Jg^Q^y obtaining 



J! 



and 



and 



where 



o(a)(y) = 3.rX„(y)ng/;(y) = ar'Al.^iy) (vf{y) \ [vf{y)M{J) 
•^oS)(y) = 9r'Al^{y)nll\y) = gr^A^^iy) {vf{y) + \ [7'J^(y) , ln(^)^ 

^tGi)(y) = Jo(Gi)(y) + iKiy) and j^YcVy) = - ^Kiy) 



|CT _ rcqp Aq jypT 

and, using Eq. H38|) . ikg can be identified as an additional, auxiliary gluon color-charge density in which 
k§(y) = -y'r'f"''r''Al,^iy)6l{y) J dzA (^^^) Al,,iz)^V^'^i.,y) . 



(E12) 



(E13) 



(E14) 



(E15) 



(E16) 



This representation enables us to express Hq — the nonlocal interaction involving quark and gluon color-charge 
densities in Eq. (|79|l — in the manifestly hermitian form 

Hc = -U dvd^dy (jo^x) + J[;|g,)(x) - *k|5(x)) D^" (r, x)92l?-(r, y) (j-(y) + Via){y) + ^Kiv) 



(E17) 

in which all operator- valued fields are hermitian as well as gauge-invariant. When we have eliminated all the 
and Hgi^^ from (-ffGi)phys and replaced them with "Pj^ and the other expressions obtained in this process, we obtain 
Eq. IHl- 



APPENDIX F 

To establish Eq. p32|l . we quote an identity in [s^ . 



ahc 



dr 



D'^"'X"'{r)Y\r)Z''{r) + X"(r)Dj''"r"(r)Z"(r) + X''{r)Y''{r)Df''Z"\r) 



0. 



(Fl) 



The proof follows from the observation that the ordinary derivative terms of the covariant derivatives, Dj 's, in Eq. (|Fip 
add up to a total derivative and the structure constant terms of Dj 's add up to zero on account of the Jacobian identity 



jlab jimc. _|_ y^^c jlma _j_ jlca jimb 



0. 



(F2) 



Notice that, the functions X , Y and Z may carry other color or vector indices and the dependence on other coordinates. 
According to Eq. (|E10(I . 



n^=(r),I?'^''(x,y) = - / dzdz'I?'^"(x,z) [n^^(x), (d ■ I?)""] <5(z - z')I?"'(z', y) 



(F3) 
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and, with Eq. ljD6|) . we have 

4(r),A^^(r)l =-[n^,^.(r),A5(r)l =-i?,Jr)[n^"(r),A5(r)" 

dxLif X>*'(x, r) dj ^'""(x, r) 9^ ^^P'^^Cx, r), (F4) 



where the symmetry property Eq. H37|l is employed to obtain the second term on the right hand side. Upon relabehng 
the dummy color indices, we have 



The 2nd term of r. h. s. of Eq. ^ = ^f^bjcmn 



dxDf I?'^-'(x,r) dj I?""'(x,r) dj a,I?"''(x,r) 



= -i^f'^bjcmbQ^-j^ic^^^ r)5jl?™(r, r), 
where the last step follows from the identities Eqs. (|F1(I and H3()() . We have, then 



(F5) 



[i;)(r), A5(r)] = _g2(2/'fc™/'^" _ /'™/"'")ajI?''''(r, r)9jl?™"(r, r) . 



Combining it with 



(F6) 



A5(r)A5(r) = -52/'''V™"9j2?''^r, r)ajl?™"(r, r) 



(F7) 



according to Eq. (|126|l and using the Jacobian identity Eq. (jF2|) . we end up with Eq. (|132|l . 
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